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Abstract. 

Context. Current models of the formation of ice giants attempt to account for the formation of Uranus and Neptune within the 
protoplanetary disk lifetime. Many of these models calculate the formation of Uranus and Neptune in a disk that may be several 
times the minimum mass solar nebula model (MMSN). Modern core accretion theories assume the formation of the ice giants 
either in situ, or between ~ 10-20 AU in the framework of the Nice model. However, at present, none of these models account 
for the spin properties of the ice giants. 

Aims. Stochastic impacts by large bodies are, at present, the usually accepted mechanisms able to account for the obliquity of 
the ice giants. We attempt to set constraints on giant impacts as the cause of Neptune's current obliquity in the framework of 
modern theories. We also use the present orbital properties of the Neptunian irregular satellites (with the exception of Triton) 
to set constraints on the scenario of giant impacts at the end of Neptune formation. 

Methods. Since stochastic collisions among embryos are assumed to occur beyond oligarchy, we model the angular momentum 
transfer to proto-Neptune and the impulse transfer to its irregular satellites by the last stochastic collision (GC) between the 
protoplanet and an oligarchic mass at the end of Neptune's formation. We assume a minimum oligarchic mass nij of 1 mg,. 
Results. From angular momentum considerations, we obtain that an oligarchic mass m, ~ Im^ < m, < 4 would be required 
at the GC to reproduce the present rotational properties of Neptune. An impact with m, > 4 is not possible, unless the impact 
parameter of the collision were very small. This result is invariant either Neptune had formed in situ or between 10-20 AU and 
does not depend on the occurrence of the GC after or during the possible migration of the planet. From impulse considerations, 
we find that an oligarchic mass m, ~ Im© < < 1.4 at the GC is required to keep or capture the present population of 
irregular satellites. If w, had been higher, the present Neptunian irregular satellites had to be formed or captured after the end 
of stochastic impacts. 

Conclusions The upper bounds on the oligarchic masses (4 m^, from the obUquity of Neptune and 1 .4 from the Neptunian 
irregular satellites) are independent of unknown parameters, such as the mass and distribution of the planetesimals, the location 
at which Uranus and Neptune were formed, the Solar Nebula initial surface mass density, and the growth regime. If stochastic 
impacts had occurred, these results should be understood as upper constraints on the oligarchic masses in the trans-Satumian 
region at the end of ice planet formation and may be used to set constraints on planetary formation scenarios. 

Key words. Planets and satellites: general - Planets and satellites: formation- Solar System: general-Solar System: formation 



1. Introduction 

The origin of the rotational properties of the planets in the Solar 
System is one of the fundamental questions in cosmogony. 
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When planetesimals are accreted, they deliver rotational angu- 
lar momentum caused by their relative motion with respect to 
the protoplanet. The angular momentum, L, acquired in an in- 
dividual collision can be in any direction. However, because 
of the symmetry of the system about the plane of the planet's 
orbit, (z, z) (-Z, -z), there is no systematic preference of 
positive or negative or L,.. An ordered component to L, is 
possible, thereby producing a net planetary spin angular mo- 
mentum either in the same direction as the orbital angular mo- 
mentimi (prograde rotation) or in the opposite direction (retro- 
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grade rotation). The angular momentum accreted from an uni- 
form and dynamically cold disk of small planetesimals results 
in a slow systematic (ordered) component of planetary rotation 
in the retrograde direction (Lissauer & Kary 1991). Lissauer 
et al. (11997b have shown that systematic prograde rotation can 
be achieved if the disk density profiles are imposed such that 
the surface mass density near the outer edges of a protoplanet's 
feeding zone is significantly greater than that in the rest of the 
accretion zone. Schlichting & Sari (i2007) obtained that plan- 
etesimals close to a meter in size, are likely to collide within 
the protoplanet's sphere of influence, thereby creating a pro- 
grade accretion disk around the protoplanet. The accretion of 
such a disk results in the formation of protoplanets spinning in 
the prograde sense. However, these models account for a net 
component of the planetary spin, and the problem of the obliq- 
uity (the angle between the rotational axis of a planet and the 
perpendicular to its orbital plane) of the ice planets remains 
open. 

Several theories have been proposed to account for the 
obliquity of giant planets. A random component of planetary 
rotation may be due to stochastic impacts with large bodies and 
may be in any direction (Lissauer & Safi'onov ll991l Chambers 
l20OT1 l. Greenberg ((T974] | and Kubo-Oka & Nakazawa (1995J 
investigated the tidal evolution of satellite orbits and examined 
the possibility that the orbital decay of a retrograde satellite 
leads to the high obliquity of Uranus, but the high mass re- 
quired for the hypothetical satellite makes this very implausi- 
ble. An asymmetric infall or torques from nearby mass concen- 
trations during the collapse of the molecular cloud core leading 
to the formation of the Solar System could twist the total angu- 
lar momentum vector of the planetary system. This twist could 
generate the obliquities of the outer planets (Ti'emaine ll991l l. 
This model has the disadvantages that the outer planets must 
form before the infall is complete and that the conditions for 
the event that would produce the twist are rather strict. The 
model itself is difficult to be quantitatively tested. In the Nice 
model of Tsiganis et al. ( I2005I I. close encounters among the gi- 
ant planets produce large orbital eccentricities and inclinations 
that were subsequently damped to the current value by gravi- 
tational interactions with planetesimals. The obliquity changes 
because of the change in the orbital inclinations. Since the incli- 
nations are damped by planetesimals interactions on timescales 
that are much shorter than the timescales for precession due to 
the torques from the Sun, especially for Uranus and Neptune, 
the obliquity returns to low values, if it is low before the en- 
counters (Hoi et al. l2007l l. Boue & Laskar (1201 Oi l reported nu- 
merical simulations showing that Uranus 's axis might be tilted 
during the giant planet instability phase described in the Nice 
model, provided that the planet has an additional satellite and a 
temporary large inclination being the satellite ejected after the 
tilt. However, the required satellite is too massive. For Saturn, a 
good case can be made for spin-orbit resonances (Hamilton & 
Ward l2004] l. but giant impacts in the late stages of the forma- 
tion of Uranus and Neptune remain the plausible explanation 
for the obliquities of the ice giants (Lee et al. l2007l l. 

If giant impacts are responsible for the obliquities of the 
ice planets, such impacts will strongly affect the orbits of hy- 
pothetical preexisting satellites. The impulse imparted at any 



collision would have produced a shift in their orbital veloc- 
ity. Satellites on orbits with too large a semimajor axis escape 
from the system (Parisi & Brunini 1997), while satellites with 
a smaller semimajor axis may be pushed to outer or inner or- 
bits acquiring greater or lower eccentricities depending on the 
impactor mass and velocity, the initial orbital elements, the ge- 
ometry of the impact and the satellite position at the moment of 
impact. The present physical and dynamical properties of irreg- 
ular satellites may be used to set constraints on the scenario of 
giant impacts at the end of the formation of ice giants (Brunini 
et al. l2052l Maris et al. l2007l Parisi et aL l2008l l. 

We investigate whether our results obtained for Uranus 
(Parisi & Brunini [T9971 Brunini et aL |2002l Parisi et aL |2008T l 
are similar when an improved method is applied to Neptune. 
Since stochastic collisions among embryos are assumed to oc- 
cur beyond oligarchy, we model the last stochastic collision 
(hereafter, GC) between the protoplanet and an oligarchic mass 
computing the angular momentum transfer to proto-Neptune 
in section 2 and the impulse transfer to its irregular satellites 
(with the exception of Triton) in section 3. In section 4, we 
describe and discuss modern scenarios of planetary formation 
in the outer Solar System. The conclusions of the results are 
presented in section 5. 

2. The spin of Neptune: Angular momentum 
transfer to Neptune by the last giant collision 
(GC) 

Beyond oligarchy, the final stage of planet formation consist of 
close encounters, collisions, and accretion events among oli- 
garchs. Strong impacts deliver spin angular momentum to the 
final planet in a random-walk fashion (Lissauer & Safronov 
119911 ). The planetary spin accumulated by successive collisions 
with a distribution of small or/and large planetesimals requires 
ad-hoc assumptions about unknown properties of the planetes- 
imal disk, such as the mass distribution of the bodies, the ve- 
locities distribution, and the regime of growth. We avoid the 
necessity of quantifying these unknown parameters modeling 
what happened to the planet just before it acquires its present 
rotational status, which is our available data. 

The last off-centre giant collision (GC) between proto- 
Neptune and the last colliding oligarch is computed assuming 
that the present spin properties of Neptune are acquired by the 
GC. The rotational parameters of the proto-Neptune prior to the 
GC are random; i.e., the rotational period and the spin obliquity 
of Neptune before the GC are unknown and taken as initial free 
parameters. The procedure by Parisi et. al. (2008) is improved 
here, since no hypothesis is applied to the spin of Neptune prior 
to the GC. We also take the velocity of the impactor as a free 
parameter that is constrained from simple dynamical bounds. 

From angular momentum conservation, we get the follow- 
ing relation between the impactor mass m, and its incident 
speed V, (Parisi & Brunini 1 19971 1. assuming that the impact is 
inelastic (Korycansky et al. 11 9901 1 : 
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where b is the impact parameter of the colHsion, £1 the present 
spin angular velocity of Neptune, £2o the spin angular velocity 
that Neptune would have today if the GC had not occurred, and 
a is the angle between Q and Qq- We take a minimum value 
of 1 ;«e for the oligarchic mass to impact proto-Neptune and 
assume a maximum oligarchic mass of 4 m@ for comparison. 

We get Neptune data from the JPL homepage. The current 
radius of Neptune is taken as 24,764 km and the mass of 
Neptune after the GC, (m,- + m^v), is taken as its current mass 
of 102.44 X 10^"^ kg. The spin period of Neptune is r= 16.11 
h, thus Q= 1.08338 x lO""* s"' (Q = In/T). Neptune's cuiTent 
obliquity is 29.58". In the single stochastic impact approach a= 
29.58°. The ice to rock ratio of Neptune is 85 %. Assuming that 
all the ice and rock are contained at the core of the planet at the 
time of the GC, the core mass of Neptune is ~ 87.074 x 10^"^ 
kg. We then compute the core radius of Neptune assuming that 
its core density and that of Uranus p^ are the same at the 
end of their formation: 
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where m„c= 74.737 x 10^^ kg is the core mass of Uranus (Parisi 
& Brunini [T997] ) and Rue = 1.8 x lO'* km its core radius 
(Bodeheimer & PoUak I1986I I. From Eq. (|2]), we get the core 
radius of Neptune, Rc = 1.9 X 10^ km. At the end of the forma- 
tion of Uranus and Neptune, their gas envelopes extend until 
the accretion radius (e.g. Bodenheimer & Pollack [T986l) . Then, 
the radius of Neptune is ~ 2.724 x 10^ km (1,100 Rn), whereas 
its core containing most of the mass has a radius of only 1.9 x 
10"* km. In this situation, a collision onto the core is necessary 
for an inelastic collision to occur and to impart the required 
angular momentum (Korycansky et al. 1990 ). Since b is an un- 
known quantity, we take its most probable value: b={2/3)Rc 
(Parisi et al. lSOUSl l. 

To set constraints on the impactor speed v,-, we calculated 
the lower bound of v, (v,„,) corresponding to a body within the 
Hill radius Rh of Neptune and undergoing a free fall towards 
Neptune's core: 
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The second term of Eq. Q is negligible, and then, v,„, depends 
only slightly on m,. Taking the average of v,,,, for m, between 
1-4 me, we get v,„, ~ 24.758 km s"'. 

The upper bound of v, (v,m) for a body bound to the Solar 
System corresponds to a parabolic orbit with the impactor ly- 
ing on the same orbital plane as proto-Neptune and moving in 
a direction opposite to Neptune's motion, including the accel- 
eration caused by the planet (Parisi & Brunini , 1997 J : 
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where Mg is the mass of the Sun and the Neptunian orbital 
semiaxis at the time of the GC. Assuming that the GC with 
Neptune occurred at 30 AU, v,m is ~ 30 km s"^ If Neptune had 



been formed between 12 and 30 AU (Dodson-Robinson et al. 
120101 Benvenuto et al. 120091 ) and the GC had occurred before 
or during outward migration (Tsiganis et al. 2005 ), v,m would 
had been between 35 and 30 km s It should be noted that 
our estimate of v,m in Eq. Q must be taken with care since 
big bodies velocity dispersion may increase so much that some 
fraction of them in the outer Solar System are ejected. Then, 
we take a maximum possible value for Vj^ of ~ 40 km s"' . 

We computed v,- as a function of Tq {To= 27r/Qo) through 
Eq. ([T]) for m, 1, 2.7, and 4 m®. Since a is a free parameter, 
we took six values of a: 0", 29.58°, 60°, 90°, 130°, and 170°. 
For a between 180° and 360°, the results would be the same as 
for the interval [0°,180°] since Eq. ([T]i is an even function of a. 
The cases with a > 60° are less probable since it requires a very 
high initial obliquity before the GC. For each a, Vi(To) must fall 
between the bounds given by v„„ and v/m- In Fig.l, we show 
the results for m,= 1 m®. The permitted values of Tq for a=0° 
(Fig.la) are between 6.95 h and 8.80 h and r,, > 43.70 h. If the 
present obliquity of Neptune was caused by a single stochastic 
impact (Fig. lb), the permitted values of To are between 7.87 h 
and 11.14 h, and Tq > 34.52 h. In Fig.lc, To must be greater 
than 12.29 h. If a > 90°, Tq must be greater than 30 h. The 
curves shift up and right as a increases. In Fig. 2, we display 
the same as in Fig.l for m,= 2.7 m^. The intervals of permitted 
To are [3.15, 4.49] h for 0° (Fig.2a) and [6.28,15.66] h for a 
- 170° (Fig.2f). For intermediate a, the range of permitted To 
falls between both cases. In Fig. 3, we show the results for m,= 
4 me, where we can see that the permitted To are very lowQ 
Then, we calculate the break-up speed given by 



- QqijRi^o, 



where /?m)=^a'/(1 + ini/3mN) (Parisi & Brunini [1 9971 1 . 

From Eq. (j5]l, we get the period Tgb (27r/Qofe) tabulated in 
Table 1. If To < Tot, the planet breaks up since centrifugal 
forces exceed gravitational forces. We then discard the cases 
shown in Fig.3a-c since the permitted Tq are very close to Tgti= 
3.43 h. It should be noted that the curves of v, in Figs. 1-3 shift 
to the left side as m, increases. Also, the cases for a > 60° are 
unlikely, since a very high initial obliquity would be required. 
Then, for an impactor mass m, > 4 me the permitted To di- 
minishes, while Toh increases (see Table 1). This implies that 
an impactor mass m, > 4 me would refute the GC hypothesis 
since Tq is less than Tob unless a is very large or b « 2/3 
Rc- However, the probability that b < Rc /3 is 0.1 and may be 
discarded. For b between Rc/3 and 2/3Rc, the probability is 
0.33, and it cannot be discarded. However, for these values of 
b, the curves of Fig. 3 shift upwards by a factor between 1 and 
2, and the permitted To remain lower than five hours for Fig. 3a- 
b. Thus, even for b between Rc/3 and 2/3Rc, the permitted Tq 



' In all the cases, v, must be zero for a= 0° and T = Tq. However, 
looking at the figures, the root of Eq. l[TJ shifts slightly to the left side 
of Tq = T as nii increases. In the development of Eq. iQ, we con- 
sidered an expression to first order in terms of the protoplanet radius 
increment for the mass increment ra,/m^ (Parisi & Brunini 1997J. For 



puting niilnifj, while for lUj- 
for m, = 4 m^. 
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2.7 rua, this error is ~ 0.02 and ~ 0.05 
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Fig. 1. Impactor incident speed v, as a function of the initial period Tq for m, - 1 m^, obtained through angular momentum 
transfer for the most probable impact parameter b = 2/3 Rc. a-a- 0", h-a- 29.58°, c-a= 60°, d-a= 90 ", e-a= 130 and f-a- 
170 The upper constraint on v,- (v,m) is depicted by dotted lines while the lower constraint on v, (v„„) is shown by dashed 
lines. v,M and v„„ are obtained from simple dynamical constraints. 
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Table 1. Break-up period of proto-Neptune before the GC for different impactor masses. 
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are close to Tot unless a is large. Moreover, if the impactor is 
large, the most probable impact parameter \s b - 2/3{Rc + R) 
(R is the impactor radius). Then, the probability that b were 
less than 2/3 Rc would be even lower. 

We may then conclude that an impact with a mass higher 
than 4 probably could not have occurred since it would be 
improbable that it could reproduce the present rotational prop- 
erties of Neptune. It might be understood as an upper bound for 
the oligarchic masses in the trans-Saturnian region at the end of 
Neptune formation. 

3. Irregular satellites: May they set constraints on 
giant impacts? 

3.1. Their origin: iiow and wiien irregular satellites of 
giant planets were captured remains open 

Rich systems of irregular satellites of the giant planets have 
been discovered very recently. The study of their origin is very 
important because it can provide constraints on the formation 
process of giant planets and may probe the properties of the 
primordial planetesimal disk (Parisi et al. 120081 Vokrouhlicky 
et al. [2008b . Neptune has at present seven irregular satellites: 
I Triton, II Nereid, IX Halimede, XI Sao, XII Laomedeia, 
XIII Neso, and X Psamathe (Holman et al. 120041 Sheppard et 



al. I2006bb . Some of their physical and orbital parameters are 
shown in Table 2. 

Irregular satellites of giant planets are characterized by ec- 
centric (highly tilted with respect to the parent planet equato- 
rial plane) and in some case retrograde, orbits. These objects 
cannot have formed by circumplanetary accretion as the regu- 
lar satellites but they are likely products of an early capture of 
primordial objects from heliocentric orbits, probably in associ- 
ation with planet formation itself (Jewitt & Sheppard 120051) . 
It is possible for an object circling the sun to be temporar- 
ily trapped by a planet. In terms of the classical three-body 
problem, this type of capture can occur when the object passes 
through the interior Lagrangian point, L2, with a very low rel- 
ative velocity. But, without any other mechanism, such a cap- 
ture is not permanent and the objects will eventually return to 
a solar orbit after several or several hundred orbital periods. 
To turn a temporary capture into a permanent one requires a 
source of orbital energy dissipation and needs for particles to 
remain inside the Hill sphere long enough for the capture to 
be effective. Although giant planets currently have no efficient 
mechanism of energy dissipation for permanent capture, sev- 
eral mechanisms may have operated at their formation epoch: 
1) gas drag in the solar nebula or in an extended, primordial 
planetary atmosphere or in a circumplanetary disk (Pollack et 
al l 19791 Cuk & Burns l2003T l. 2) Pull-down capture caused by 
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the mass growth and/or orbital expansion of the planet, which 
expands its Hill sphere (Brunini [l995l Heppenheimer & Porco 
I1977I I. 3) Collisional interaction between two planetesimals 
passing near the planet or between a planetesimal and a reg- 
ular satellite. This mechanism, the so-called break-up process, 
leads to the formation of dynamical groupings (e.g. Colombo 
& Franklin 1971, Nesvorny et al. I2004I I. After a break-up the 
resulting fragments of each progenitor form a population of ir- 
regular satellites with similar surface composition, i.e. similar 
colors and iiTegular shapes, i.e. light curves of wide amplitude 
(Maris et al. 120011 Maris et al. 120071) . 4) Collisionless inter- 
actions between a massive planetary satellite and guest bodies 
(Tsui |1999l l or between the planet and a binary object (Agnor & 
Hamilton 120061 Vokrouhlicky et al. l2008l l. 5) Scattering by an 
eccentric Triton applicable to the Neptunian system (Goldreich 
et al l 19891 Cuk & Gladman I2005I I. and 6) collisionless inter- 
actions within the framework of the Nice model, where capture 
occurs during migration through three-body interactions dur- 
ing close encounters between the giant planets (Nesvorny et 
al. |2007, Bottke et al. I2010I I. There might have been several 
encounters and stochastic collisions among protoplanets dur- 
ing migration, all of which could have both removed and cap- 
tured iiTegular satellites. It was shown that the last captured 
population of irregular satellites by three-body encounters in 
the framework of the Nice model were very likely significantly 
larger than what we observe today. These populations had then 
experienced rapid collisional evolution and almost self destruc- 
ted, evolving into quasi-steady state size frequency distribu- 
tions with extremely shallow power-law slopes for D> 10 km 
(Bottke et al. THW). 

The knowledge of colors and of the size and shape distri- 
bution of iiTegular satellites is important for knowing their re- 
lation to the precursor population. It brings valuable clues to 
investigate if they are collisional fragments from break-up pro- 
cesses occurring at the planetesimal disk and thus has nothing 
to do with how they were individually captured later by the 
planet, or if they are collisional fragments produced during or 
after the capture event (Nesvorny et al. 120031 12007i Bottke et 
al. 120 1 Oi l. Despite many studies in past years, inconsistencies 
are present between colors derived by different authors for the 
in-egular satellites of Uranus and Neptune (Maris et al. 120011 
Maris et al. 120071 Grav et al. I2004I I. Moreover, if the precur- 
sors of the Uranian and Neptunian irregular satellites are ob- 
jects of the Kuiper Belt or C-,P-, and D-type asteroids is still a 
matter of debate (Bottke et al. |20I0l Levison et al. ,2009) . Even 
the recent proposed capture models run into trouble using the 
size frequency distribution of irregular satellites as constraints 
(Nesvorny et al. 2007, Bottke et al. ,20 10 J . An intensive search 
for fainter irregular satellites and a long term program of ob- 
servations to recover light curves, colors, and phase-effect in- 
formation in a self-consistent manner is needed. Constraints on 
the different dissipative mechanisms for the permanent capture 
of iiTegular satellites should be investigated in the framework 
of the different models proposed for the formation of the ice 
giants. 

Although, recent models study the capture of irregular 
satellites by close encounters between giant planets in the 
framework of the Nice model, there is very httle work on the 



capture and loss of irregular satellites by stochastic collisions 
(Parisi et al. I2008I I. Here, we investigate whether the present 
population of Neptunian irregular satellites could have been 
captured by stochastic impacts and/or if they could have sur- 
vived to stochastic impacts. We attempt to set constraints on 
the occurrence of giant impacts at the end of the formation of 
the ice giants, on the impactor masses, and on the epoch of the 
capture/origin of the Neptunian irregular satellites. 

3.2. The Neptunian irregular satellites were captured 
before, after, or during stochastic impacts? 

All stochastic impacts add angular momentum to the random 
component of planetary spin and transfer impulse to the planet 
and its satellite system. Although the transfer of impulse at any 
collision onto the planet affects the whole system, the trans- 
fer of angular momentum affects only the rotational properties 
of the planet but not the orbital or rotational properties of its 
satellites. 

The impulse accumulated by successive collisions with a 
distribution of small or/and large planetesimals requires ad- 
hoc assumptions about the unknown planetesimal mass distri- 
bution, the unknown initial surface mass density, and the un- 
known regime of growth. In the GC scenario, we avoid the ne- 
cessity of quantifying these parameters, since we assume that 
the orbital properties of the irregular satellites of Neptune be- 
fore the GC are unknown as are the physical and dynamical 
properties of oligarchs. These initial free parameters are con- 
straint by the model from the current orbital properties of the 
irregular satellites and using simple general dynamical bounds. 

The transfers of six of the seven irregular satellites of 
Neptune to their current orbits by the GC are computed fol- 
lowing a somewhat more improved procedure than in Parisi et 
al. ( I2OO8I 1. Irregular satellites may be pushed to outer or inner 
orbits. Moreover, they may be captured or unbound by the GC. 
We attempt to set constraints on their initial orbital properties, 
the impactor mass, and the epoch of their origin. We do not 
include Triton in our treatment, since this satellite is assumed 
to have a different origin than the other irregular satellites of 
Neptune owing to its high mass, small orbital semiaxis and zero 
orbital eccentricity (see Table 2). 

Just before the GC, the square of the orbital velocity vi of 
a preexisting satellite of negligible mass is given by 

= Gm^J^ - Ij, (6) 

where r is the position of the satellite on its orbit at the moment 
of the GC and oi its orbital semiaxis. 

After the GC, the satellite is transferred to another orbit 
with semiaxis 02 acquiring the following square of the velocity: 

vj = G(mN + mi)i- - —\ (7) 
\r aoj 

We set Vj = A and v\ - B (\ -¥ milmj^) v^, where A and 
B are arbitrary coefficients (0 < A < 1, B > 0), being the 
escape velocity at r before the GC. 
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Table 2. Orbital and physical parameters of the Neptunian irregulars taken from JPL (http://ssd.jpl.nasa.gov i. 
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Irregular Satellites 


K 


pT 


Magnitude 


Albedo 


a" 


e 


i+ 


P++ 


I Triton 


1,352.6 


2.06 


13.5 V„ 


0.76 


14.33 


0.0000 


156.885 


5.88 


n Nereid 


170.0 


1.50 


19.7 V„ 


0.16 


222.65 


0.7507 


7.230 


360.13 


IX Halimede 


31.0 


1.50 


19.2 R 


0.04 


670.77 


0.2646 


112.712 


1879.08 


XI Sao 


22.0 


1.50 


24.5 R 


0.04 


897.60 


0.1365 


53.483 


2912.72 


XII Laomedeia 


21.0 


1.50 


25.5 R 


0.04 


951.66 


0.3969 


37.874 


3171.33 


XIII Neso 


30.0 


1.50 


24.6 R 


0.04 


1990.19 


0.5714 


136.439 


9740.73 


X Psamathe 


20.0 


1.50 


25.5 R 


0.04 


1942.17 


0.3809 


126.312 


9074.30 



* The mean satellite radius r, is given in km. 

** The mean satellite density p, is given in g cm"^. 

X The orbital semiaxis a is in units of Rt^= 24,764 km. 

+ The orbital inclination i is measured with respect to the ecliptic. 

++ The orbital period P is given in days. 



The semiaxis of the satellite orbit before (ai) and after (02) 
the GC verify the following simple relations; 



2(1 -A) 



2(1-B) 



(8) 



If A < B then ai < 02- In the special case of B = 1, the orbits 
are unbound from the system. If A > Z? then a; > 02, the initial 
orbit is transferred to an inner orbit, providing a mechanism 
for the permanent capture of the iiTegular satellite even from 
an unbound orbit (A > 1). When A - B, the orbital semiaxis 
remains unchanged (ai = 02). 

The position r of the satellite on its orbit at the epoch of the 
impact may be expressed in the following form (Bunini et al. 
laofel Parisi et al. [2008) : 



2 G niN 

(Ay)2 



B' -A 



n2 



VAcos^+ yliB^'^^A)TAco^ 



(9) 



with B' - B {I + miliiiM), and 'I' is the angle between vi and 
the orbital velocity change AV imparted to Neptune by the GC. 
The value of AV is obtained through impulse conservation con- 
siderations at collision (Paiisi & Brunini ,1997j : 



(rrii + OTAr)Ay — rriiVi. 



(10) 



Upper bounds on oi (qim) and 02 («2m) are obtained from 
Eqs. (j8]l and (|9]l with *I'=180°, i.e., assuming the impact in 
the direction opposite to the orbital motion of the satellites and 
taking the positive sign of the square root: 



a\M — 



G niN (B' - AY 



(Ay)2 (1 -A)( Vb' - Va)2 

GniNiB' -Af 

(Ay)2 (1 -B)( Vb'- Va)2 



(11) 



The minimum eccentricity of the orbits before the collision 
is given by 



ei„„„=2(l-A)-l // A<0.5 
ei,„,„ = 1-2(1-A) // A>0.5, 



(12) 



while the minimum eccentricity of the orbits after the collision 
is 



e2„H« = 2(1 -B)- 1 // B<0.5 
e2m,« = 1 -2(1 -B) // B>0.5. 



(13) 



Using the current orbital elements of the Neptunian irreg- 
ular satellites as initial conditions (see Table 2), their orbital 
evolution over a period of 10^ yrs is computed using the inte- 
grator evorb 13 of Fernandez et al. (120021) . where the perturba- 
tions of the Sun, Jupiter, Saturn and Uranus are included. The 
minimum, mean, and maximum eccentricities and semiaxes are 
shown in Table |3]for all the Neptunian irregular satellites with 
the exception of Triton. We assume that gravitational perturba- 
tions are the only effect capable of altering the orbital elements 
of the iiTegular satellites during Solar System evolution after 
the GC since gas drag caused by an ice giant's envelope proves 
negUgible (Parisi et al. l2008b . 

We compute a\M and a2M from the upper and lower bounds 
of Ay obtained through Eq. ( 10 1 using v,=v,m and Vi-Vim, re- 
spectively. We take vim- 40 km s"' and v,,,, from Eq. We 
take nii- 1 ;?% as the minimum impactor mass. For each A, we 
then calculate the values of B (B = B7(l + niilm;^)) corre- 
sponding to the transfers to a2M- {cimax, o.,„in\, where a^ax and 
a,„i„ are the current maximum and minimum semiaxes of the 
irregular satellites tabulated in Table 3. Introducing these val- 
ues of B in Eq. ( fTSj l, we get the minimum possible values of 
e2mm, 62,,,, that the orbit of each irregular satellite may acquire 
at impact for each initial condition A. 

The transfer of each irregular satellite from its original or- 
bit to the present one is only possible for those values of (A, B) 
which satisfy the condition e2m < e„,a.x (emax tabulated in Table 
3). Transfers with A < B, imply that the irregular satellite be- 
longs from an inner orbit, while transfers with A > B imply 
that the iiTegular satellite belongs from an outer orbit; i.e., the 
GC may produce the permanent capture of a body undergoing 
a temporary capture or even in heliocentric orbit (A > 1). 

The transfers of Halimede and Nereid are shown in Fig.4, 
those of Sao and Laomedeia in Fig. 5, and the transfers of Neso 
and Psamathe are displayed in Fig. 6. The lower curve in Figs.4- 
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Table 3. Variation in the orbital eccentricity and semiaxis (in units of R^- 24,764 km) of the Neptunian irregulars due to Solar 



t perturbations 


over a 


period of 


10^ yrs. 




- a,ni„ and A„ 




^mean • 




Satellite 






^max 






^^max 


K ■ 




Nereid 


0.7426 


0.7496 


0.7566 




222.659 


222.687 


0.026 


0.028 


Halimede 


0.2073 


0.4331 


0.9026 


668.865 


670.297 


672.025 


1.432 


1.728 


Sao 


0.0628 


0.2839 


0.6619 


892.175 


896.893 


902.070 


4.718 


5.177 


Laomedeia 


0.2648 


0.4015 


0.5482 


942.266 


948.295 


955.526 


6.029 


7.231 


Neso 


0.0942 


0.4140 


0.8803 


1,919.790 


2,013.260 


2,130.300 


93.470 


117.040 


Psamathe 


0.0482 


0.3970 


0.9144 


1,832.740 


1,912.260 


2,015.540 


79.520 


103.280 



Halimede 



em 



Nereid 



©In 



1 O.B 0,6 0,4 0,2 0,0 0,2 0,4 0,6 0,8 1 1 0,8 0,6 0,4 0,2 0,0 0,2 0,4 0,6 0,B 1 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



0,1 0,2 0,3 0,4 0,5 OB 0,7 0,6 0,9 



1 0,2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Fig. 4. The transfers capable of producing the present orbits of Halimede (left panel) and Nereid (right panel). Lower curve 
corresponds to the minimum impactor speed v„„ and an impactor mass m, = 1 mg. Intermediate curve corresponds to the maximum 
impactor speed y,M and nij- 1 me- Upper curve corresponds to v,„, and m,= 3.5 me. The maximum current eccentricity e^ax is 
shown by a dotted line. For Halimede the upper curve and e,„a j overlap. A (B) is the square of the ratio of the satellite's speed just 
before (after) the impact on the escape velocity at the satellite's location just before (after) the impact. ei,„ (eim) is the minimum 
eccentricity of the orbits before (after) collision. 



6 represents the transfers for v,=v,>„ and m,= 1 m®. It should 
be noted that by minimizing v, and m,, we are maximizing 
the number of transfers. The real eccentricity that the irregu- 
lar satellite may acquire for each v, and m, may take any value 
between e2m and e,„ax- The probability of the existence of the 
irregular satellite prior to the GC diminishes as e2m approaches 

In Fig.4, the intermediate curve represents the transfers for 
mi= 1 me and v,=v,m- The transfers for m,- 1 and v„„ < 
V, < v,M fall between the lower and intermediate curves. The 
upper curve corresponds to v„„ and m,= 3.5 me, and e,„„j and 
the upper curve overlap. Halimede and Nereid then set an upper 
bound on m, of 3.5 m^. Within this mass limit, the GC itself 
could provide a mechanism for the capture of Halimede and 
Nereid. 

In Fig. 5, the transfers for Sao are shown in the left panel. 
There are no transfers for v,m for any impactor mass. The upper 
curve shows the transfers for v,>„ and m,= 1.6 me, which is very 
close to e,„ax ■ There is no transfer for m, > 1.7 m^. Sao set an 
upper bound of 1.7 me on the impactor mass. The transfers of 



Laomedeia are displayed in the right panel. The upper curve 
displays the transfers for and m,= 1.35 me, which are the 
same as for v/^ and m,= 1 me. Laomedeia set an upper bound 
on m, of 1.4 m^. Within this mass limit, the GC itself could 
provide a mechanism for the capture of Sao and Laomedeia. 

The transfers of Psamathe are shown in the left panel and 
those of Neso in the right panel of Fig. 6. The intermediate 
curve corresponds to the transfers for v,m and m,= 1 m^. The 
upper curve coiTesponds to v,>„ and m,= 2.3 me (Psamathe), and 
m,= 1.8 me (Neso). Psamathe set an upper bound on m,- of 2.3 
me, while Neso of 1.8 me. Within this mass hmit, the GC itself 
could provide a mechanism for the capture of Psamathe and 
Neso. The transfers of Fig.6 have been computed over a range 
of 200 Rn around the current orbital semiaxes of both irregular 
satellites (see Table 3). However, e2m lies close to e,„ax for v; > 
v,„, and m; > 1 me. This result makes the existence of Neso and 
Psamathe piior to the GC, as well as their capture by the GC, a 
possible but not very probable event. 

Parisi et al. (.2008 ) carried out a somewhat similar model of 
the Uranian system. They obtained strong constraints for two 
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Fig. 5. The transfers capable of producing the present orbits of Sao (left panel) and Laomedeia (right panel). Lower curve 
corresponds to the minimum impactor speed v,„, and an impactor mass m,= 1 m®. Upper curve left panel corresponds to v,„, and 
m,= 1.6 m^j which is the same as for v,m and m,= 1 m^. Upper curve right panel corresponds to v,„, and m,- 1.35 m^g. There is no 
transfer for v,m and m,= 1 m®. The maximum current eccentricity e,„ax is shown by a dotted line. A (B) is the square of the ratio 
of the satellite's speed just before (after) the impact on the escape velocity at the satellite's location just before (after) the impact, 
eim (eim) is the minimum eccentricity of the orbits before (after) collision. 



Uranian irregular satellites: Prospero and Trinculo. The only 
transfers for both irregular satellites were close to the pericenter 
of an eccentric initial orbit (ei,„ > 0.58 for Prospero and ei,„ > 
0.62 for Trinculo). 62,,, for Trinculo resulted in the range [0.16- 
0.23], very close to emax(0.237). Moreover, for Prospero e2m 
was ~ e„ax (0.571). This result implied that Prospero could not 
survive stochastic impacts. Either Prospero was created after 
the epoch of giant impacts or giant impacts did not occur 

In this paper, Laomedeia set an upper constraint on the im- 
pactor mass of 1.4 m^. Assuming that the Neptunian irregu- 
lar satellites belong to a common origin, then they set an up- 
per constraint of 1.4 on the impactor masses in the trans- 
Saturnian region at the end of Neptune formation. 

Assuming that giant impacts occurred beyond oligarchy, 
our results imply that either the oligarchic masses in the trasn- 
Saturnian region at the end of Neptune formation were less that 
1.4 me, or the current Neptunian iiTegular satellites had to be 
captured after the end of stochastic impacts. Alternatively, par- 
ent objects with orbits different to the iiTegular satellites' cur- 
rent orbits could have been captured by any process before or 
during stochastic impacts without restrictions on the impactor 
mass, but the present population of irregular satellites should 
then be the result of their later coUisional evolution. This pos- 
sibility would agree with the results of Bottke et al. (1201 Oi l. 

4. Discussion 

It has long been known that dynamical times in the trans- 
Saturnian MMSN are so long that core growth takes more than 
15 Myr Observations of young, solar-type stars suggest that 
circumstellar disks dissipate on a timescale of few Myr (e.g. 



Briceno et al. l200H . Different models and scenarios have been 
proposed to account for the formation of the ice giants within 
the protoplanetary disk lifetime. Modern models shorten the 
timescale for giant planet formation if taking a higher initial 
surface density into account well above that of the MMSN 
and/or the formation of all giant planets in an inner compact 
configuration (Dodson-Robinson et al. 2010 Benvenuto et al. 
120091 Thommes et aL l2003l Tsiganis et al EOOST l. 

Core accretion models of giant planet formation in situ 
(Pollack et al l 1996b and core accretion models with the gi- 
ant planets forming at the initial semiaxes of the Nice model 
(Dodson-Robinson et al,2010 Benvenuto et al.|2009jl compute 
the solid component because of the accretion of planetesimals. 
However, no core accretion model has at present, considered 
the angular momentum transfer during accretion and the pos- 
sible accretion of other protoplanets onto the core with the ex- 
ception of Korycansky et al. ( I1990I I who carried out a hydro- 
dynamical model of a giant impact on Uranus. 

Modern scenarios have several difficulties to overcome, and 
inconsistencies among the different models and scenarios are 
still present. We discuss four important issues below. 

1) Most core accretion models (e.g. Dodson-Robinson et 
al. 1201 Oi l, assume planetesimals encountering the core at the 
Hill velocity; i.e., they do not take the stirring of these bodies 
into account. The Hill velocity vh is a characteristic velocity 
associated with the Hill radius Rh and is defined as vh - ^^Rh, 
where Q is the orbital frequency around the Sun. A speed of 
the surrounding planetesimals higher than vjf would diminish 
the rate of growth of the ice giants in the massive planetesimals 
disk: either they had formed in situ or at 10-20 AU. 
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CQ 
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ViM,mi=1m0 ^im:mi=1.8m0 
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1 

■ 0,8 

■ 6 

■ 0,4 

■ 02 

■ 0,0 

■ 0,2 

■ 04 
■0,6 
,0,8 



Fig. 6. The transfers capable of producing the present orbits of Psamathe (left panel) and Neso (right panel). Lower curve cor- 
responds to the minimum impactor speed v„„ and an impactor mass m,= 1 m^. Intermediate curve corresponds to the maximum 
impactor speed Vj^ and m,= 1 m^,. Upper curve left panel corresponds to v,„, and m,= 2.3 m^. Upper curve right panel corresponds 
to Vim and m,= 1.8 me. The maximum current eccentricity e,„av is shown by a dotted line. In both panels the upper curve and 
Smax overlap. A (B) is the square of the ratio of the satellite's speed just before (after) the impact on the escape velocity at the 
satellite's location just before (after) the impact, ej,,, (e2m) is the minimum eccentricity of the orbits before (after) collision. 



2) An initial mass 5-10 the MMSN is required to shorten 
the timescale of ice giants formation, but within a 10 MMSN, 
Jupiter falls like a stone into the Sun due to type III migration 
(Crida l2009l l. 

3) If oligarchic masses can reach the isolation mass in 1- 
10 Myrs at 10-30 AU is unknown. The final mass of the oli- 
garchs in the outer solar System remains an open question. 
Without collisions among oligarchs, the mass of the core of 
proto-Uranus and proto-Neptune in the frame of the MMSN is 
the mass of an oligarch; i.e., much low than the present solid 
cores of Uranus and Neptune (either the ice planets had formed 
in situ or between 10-20 AU). Collisions among oligarchs takes 
a very long time to form the ice giants cores. Then, an ini- 
tial surface density ~ 5-10 times that of the MMSN would be 
required to produce oligarchs with masses similar to those of 
the ice giants cores. However, a solid surface density of this 
size would lead to the formation of about five ice giants in- 
stead of two, which occurred with the other three giants; i.e., 
whether they were ejected or if they simply were spread out 
being all retained is a matter of debate (Goldreich et al. 2004 
Dodson-Robinson et al. 120101 Ford & Chiang 120071 Levison 
& Morbidelli 12007b . In the last case, then, where are they? 
However, Thommes et al. ( I2003I I have shown that even in a disk 
10 times the MMSN, oligarchs do not have time to reach their 
isolation mass in the outer Solar System, and even an Earth 
mass at the orbit of Uranus by 10 Myrs is implausible. 

4) Actually, the mass of the planetesimals from which 
Uranus and Neptune accreted is a matter of debate. Goldrecih 
et al. ( 1200 41) suggest that, in the dispersion dominate regime the 
velocity dispersion of small bodies u could exceed their sur- 



face escape velocity, and thus collisions at these speeds would 
be destructive. They speculate that before the oligarchs be- 
come planets, small bodies could fragment, reducing their size. 
Smaller size would imply smaller u, hence faster accretion, thus 
returning to the case u < vh &t isolation. However, there is no 
reason planetesimals should fragment in the dynamically cold 
disk required to produce Uranus and Neptune in situ (Dodson- 
Robinson et al. l2010l l. It is usually believed that the planetesi- 
mal distribution follows a power law of the type dn cc m^^dm 
(e.g. Benvenuto et al. 120091 and references therein). Indeed, 
Jupiter family comets follow a distribution of this type with 
a - 11/6. Benvenuto et al. (2009) computed the formation of 
Uranus and Neptune assuming the initial semiaxes distribution 
of the Nice model taking a - 2.5. They approached the con- 
tinuous distribution to a discrete distribution with planetesimal 
sizes in the range 30 m 100 km. However, Dodson-Robinson 
et al. (120101) computed the accretion of Uranus and Neptune in 
the frame of the Nice model, too, but considered planetesimals 
with radii of several hundred kilometers since they assumed the 
theory of planetesimal formation based on the streaming insta- 
bility which produces planetesimals of about 600 km (Johansen 
et al. 2007 ). Although the streaming instability has several ben- 
efits, observational evidence (e.g., comet outbursts and chon- 
drules) shows a fractal structure of the planetesimals consistent 
with coagulation processes for planetesimal formation (Paszun 
& Dominik |2009l l. According to Morbidelh et al. ( |20D9l l, as- 
teroids were born big, suggesting that the minimal size of the 
planetesimals was ~ 100 km. 

In summary, it is not known whether the ice planets formed 
in situ, or well inside the 20 UA and/or if the initial mass of 
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the nebula was that of the MMSN or much larger. Moreover, 
the mass of the planetesimals from which Uranus and Neptune 
accreted remains a matter of debate. It is necessary to look for 
independent ways of setting constraints on models of ice giant 
planet formation. 

Our models are independent of unknown parameters, such 
as the mass and distribution of the planetesimals, the location 
at which Uranus and Neptune were formed, the Solar Nebula 
initial surface mass density, and the regime of growth. Our con- 
straints on the oligarchic masses may be used to set constraints 
on planetary formation scenarios. 

5. Conclusions 

We have modeled the angular momentum transfer to proto- 
Neptune and the impulse transfer to its irregular satellites by 
the last stochastic collision (GC) between the protoplanet and 
an oligarchic mass. Since stochastic impacts are aleatory and 
their number is uncertain, the rotational properties of proto- 
Neptune and the orbital properties of the Neptunian system be- 
fore the GC are taken as initial free parameters in our model 
and are constrainted using the present rotational parameters 
of Neptune and the present orbital and physical properties of 
Neptune and its irregular satellite population. 

Assuming a minimum impactor mass m,~ 1 »%, the mass 
of the last impactor m, < 4 77% is required to account for the 
current spin properties of Neptune. This result is invariant: 
Neptune had formed either in situ or between 10-20 AU and 
does not depend on the possible occuiTence of the stochastic 
impact during or after migration. A collision with m, > 4 
could not have occuiTed since it cannot reproduce the present 
rotational and physical properties of the planet, unless the im- 
pact parameter of the collision was very small. The formation 
of Neptune as the result, for instance, of coUisional accretion 
between two similar oligarchs with masses ~ 7 seems to 
be unlikely. Assuming the occurrence of giant impacts as the 
cause of Neptune's obliquity, the 4 m® mass limit must be un- 
derstood as an upper bound for the oligarchic masses in the 
trans-Saturnian region at the end of Neptune's formation. 

The study of the origin of irregular satellites of giant planets 
is very important because it puts constraints on formation pro- 
cesses of giant planets and may probe the properties of the pri- 
mordial planetesimal disk from which irregular satellites were 
captured. Although recent models study the capture of iiTegu- 
lar satellites by close encounters between giant planets in the 
framework of the Nice model (Bottke et al. 1201 Oi l, there is very 
little work on their capture and loss by stochastic collisions 
(Parisi et al. l2008T l. 

We found that either the mass of the last impactor on 
Neptune was less than ~ 1 .4 m^, or the present Neptunian ir- 
regular satellites had to be formed or captured after the end of 
stochastic impacts. This result is invariant: Neptune had formed 
either in situ or between 10-20 AU and does not depend on the 
possible occurrence of the stochastic impacts during or after 
migration. 

If the present population of irregular satellites were cap- 
tured after the end of stochastic impacts, the mechanism of 
capture able to operate at such late stages in the scenario of the 



formation of the planet should be investigated. Alternatively, 
parent objects with orbits different from the irregular satellites' 
current orbits could have been captured by any process be- 
fore or during stochastic impacts, but the present population 
of irregular satellites should be the result of their later coUi- 
sional evolution. This possibility would agree with the results 
of Bottke et al. (| 20T0] |. 

Colors are an important diagnostic tool in attempting to un- 
veil the physical status and the origin of the Neptunian irregular 
satellites. In particular it is very important to assess whether it 
is possible to define subclasses of irregular satellites just by 
looking at colors and by comparing colors of these bodies with 
colors of minor bodies in the outer Solar System. The knowl- 
edge of the primordial population of iiTegular satellites could 
offer valuable clues to know whether Neptune was formed in 
situ or between 10-20 AU. An intensive search for fainter ir- 
regular satellites and a long term program of observations able 
to recover lightcurves, colors, and phase effects in a "self con- 
sistent" manner is mandatory. 
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